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We present analytic expressions that approximate the behavior of the spacetime of 
a collapsing spherically symmetric scalar field in the critical regime first discovered by 
Choptuik. We find that the critical region of spacetime can usefully be divided into a 
\^ , "quiescent" region and an "oscillatory" region, and a moving "transition edge" that 

^\ • separates the two regions. We find that in each region the critical solution can be 

^\ ' well approximated by a fiat spacetime scalar field solution. A qualitative nonlinear 

, matching of the solutions across the edge yields the right order of magnitude for the 

f— I ' oscillations of the discretely self-similar critical solution found by Choptuik. 

: 

o\ : 

T— I I. INTRODUCTION 

> ; 

■ The dynamics of a spherically-symmetric massless scalar field coupled to general relativity can have 
, two late-time outcomes. Either a black hole forms or the scalar field energy eventually radiates away, 
T— I ' leaving spacetime flat. By studying solutions very close to the critical case (formation of a black hole 
^) ' of vanishing mass) Choptuik has found two very interesting phenomena, as yet not explained in 
any fundamental way: universality and echoing. Similarities to the Choptuik results have been found 
in other contexts. Abrahams and Evans j^] studied the collapse of axisymmetric vacuum gravitational 
O • fields. Coleman and Evans pi have studied the collapse of spherically symmetric radiative fluids, 
O^" and Eardley and Hirshmann [W studied the collapse of spherically symmetric complex scalar fields. 
^ ', "Universality" is a feature of supercritical cases, those in which holes form. The dependence of 
the mass of the final hole on a parameter of the initial data, appears to have a universal power-law 
^ character and this character, as well as the numerical value of the exponent appear to be the same 
for different families of data, though not, as originally thought, for different collapsing fields 
' This feature of gravitational collapse would seem, therefore, to have the same universal character 
— as well as the same mathematical appearance — as phase transitions in statistical mechanics. 

In sharp contrast, "echoing" defies close comparison with any other physical phenomenon. For 
solutions very close to criticality, the evolving fields go through a sequence of cycles, each of which 
duplicates the previous one, but on a radial scale smaller by around a factor of 30, and on a time scale 
smaller by the same factor. This has been seen only in the work of Choptuik 0| and Abrahams and 
Evans |^ axially symmetric examples, in which the echoing occurred with a scaling factor around 6, 
rather than 30 as in Choptuik's spherically symmetric results. Very recently Gundlach |Q, with an 
independent numerical approach, has confirmed that the echoing is an exact feature of the solution. 
For fluid collapse |^ and the complex scalar fleld Q exact self-similar critical solutions arc known, 
with no echoing (in the complex scalar fleld case there is echoing in the phase) . It is widely suspected 
that the echoing behavior has an origin related to self-similarity, but a detailed understanding is still 
lacking. 

In this paper we will consider only that feature of Choptuik's solution which we consider most 
unusual, and surely that which is least understood: the echoing. Since no "simple" model problem 
is yet know which has this feature, we focus instead on the equations of Choptuik's physical problem 
to see the manner in which the echoing solution arises. We are, then, in some sense dealing only 
with the "how" of the solution, rather than the "why," in the hope that this might be a useful step. 

With this aim we will divide the spacetime into two regions where we are able to find approximate 
solutions to the coupled Einstein-Scalar field equations and we will discuss a nonlinear matching 
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between them. It wil turn out, somewhat surprisingly, that in the two relevant regions, the behavior 
is mainly dominated by flat spacetime evolutions. The oscillatory behavior, which ab initio appears 
as a feature of the nonlinearities of general relativity, can actually be accounted for by flat spacetime 
behavior. The nonlinear nature of the theory is necessary to join the solution in the two regions, 
and this join determines the scaling factor of the echoing. 

The plan of this paper is as follows. In section II we present the basic equations and discuss the 
results for the dynamics observed by Choptuik. In section III we present our approximate solutions 
and discuss their regions of validity. In section IV we discuss the relevance of these approximate 
solutions. 



II. BASIC EQUATIONS AND CHOPTUIK'S SOLUTION 

The physical problem is the dynamics of a scalar field (jj coupled to general relativity. The combined 
set of Einstein-Scalar equations are, 

R^ty = 87r9^(/)9^0 (1) 
□0 = 0. (2) 

We assume spherical symmetry for the scalar field and the spacetime, so the metric can be written 
as 

ds^ = ^a"^ {r,t)dt^ +a^{r,t)dr^ +r^{d0^ +sm'^e d(p'^) . (3) 
It is useful to introduce the mass function 

m{r,t)=r (4) 

and the dimensionless "phase space" variables used by Choptuik: 

X(r,t) = V2^^^, Y{r,t) = V2^'-^. (5) 
a or am 

With these choices, the Einstein-scalar equations reduce to a set of four partial differential equations, 

{aX),t^-— + {aY).r (6) 
r 

(ar),, ^ — + {aX)^r (7) 
r 

(7)^ + 5 = J(^' + ^') (8) 
aa r 

The first two equations are equivalent to the Klein-Gordon equation in curved spacetime, and the 
last two are the Einstein equations determining the metric coefficients a and a as functions of the 
energy density of the scalar field, which is proportional to + Y^. 

Note that the form of these equations is invariant under the transformation 

a{r,t) a{r,t)a{t) dt ^ dt/a{t) . (10) 

This corresponds to the fact that the t coordinate is not completely fixed by the form of (||). The 
degree of freedom inherent in a{t) can be eliminated by fixing the function a{r, t) at one value of r. 
In this paper we make the choice a{r = 0, i) = 1 for all t, so that the t coordinate measures proper 
time of an observer at the origin. 

For near critical solutions, there is a sequence of echoing oscillations, each shorter than the previous 
one by a scaling factor of around 30. For supercritical solutions the sequence ends with horizon 
formation; for subcritical solutions the echo sequence dies out. For a perfectly critical solution, the 
solution studied by Gundlach |^], the echoing would go through an infinite number of cycles. In all 
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cases, the duration of proper time is finite, and is not highly sensitive to the way the sequence ends. 
From computational results with many echoes, then, the time for the end of the sequence can be 
estimated accurately. We adjust our time coordinate, so that t» = 0. That is, our zero of t is chosen 
to be the endpoint of echoing. 

Having fixed the zero and the scale of our t coordinate, we next define the logarithmic coordinates 



Inr 



H-t) 



(11) 



Notice that our choice of a, and hence of t, differs from that of Choptuik ^ [he chooses 
a{r = oo,t) = 1] but p, r defined in ( [III ) agree with Choptuik's logarithmic coordinates p, r. 

Figure 1 shows a series of snapshots of the scalar field variable X as a function of p taken at 
various late times r. These are data made available by Choptuik on the internet from his numerical 
evolutions. In each of the plots we can distinguish two regions. On the left of the plots there is a 
"quiescent" region in which the solution dies off; on the right is an "oscillatory" region in which the 
solution has a sinusoidal appearance. We call the moving boundary between the two regions the 
"transition edge." 
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FIG. 1. The Choptuik numerical results for the scalar field variable X as a function of the radial 
coordinate p, at different times r. 



"Echoing" of the solution is the property that the pattern of X{r,t) repeats itself on a changed 
radial and time scale according to 



or equivalently 



X{r,t) =X{e\,e^t), 
X{p,T)=X{p + A,r + A) 



(12) 
(13) 
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If A were a continuous variable, this property would tell us that X is "self similar," i.e., a function 
only of r/t. But equations (^2|) and ( |l3[) are valid only for a fixed numerical value for A, approxi- 
mately 3.4 (a value of 3.4439 ± 0.0004 according to Ref. Q). For this reason the echoing property 
has been referred to as a "discrete self similarity." 

Figure 2 shows the behavior of other variables of the problem, as functions of p at a fixed r. 
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FIG. 2. The Choptuik numerical results for scalar field and metric variables, as a function of p, all at 
time T = -7.87488. 

The scalar field variables X and Y have the appearance of being sinusoids approximately 150 
degrees out of phase. Choptuik ||^ has shown that it is useful to view the relation of X and Y in 
phase space plots, as in Fig. 3. 

The elliptical portion of the plot corresponds to the oscillatory region of the solution. The tran- 
sition region appears as the vertical spike which connects the elliptical windings to the origin. As 
radius increases, the windings decrease in size and the curve connects to the origin at a point rep- 
resenting spatial infinity in the numerical solutions. In a truly critical solution there would be an 
infinite number of windings and the connection with spatial infinity would be absent. 
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FIG. 3. The Choptuik numerical results for the relationship of the phase space variables X and Y . Each 
plot corresponds to a single moment of t and the full range of radii. The elliptical winding spiralling into the 
origin corresponds to radius growing large; the vertical spike approaching the origin corresponds to radius 
going to zero. 

The qualitative features of a, in Fig. 2, can be understood from (|^). In this equation a^p is "fed" 
by the source term a^{X^ + Y"^) which at each instant of time is a quantity bounded between zero 
and a value less than one. To get a rough understanding of the behavior of a we can take the source 
term as approximately constant and then the solution for a is 

a - T{T)e''P , (14) 

where fc is a constant related to the average value of a^(X^ + Y"^) and is a function of time. 
In view of the echoing, it is worthwhile to rewrite the equations in term of a new variable /i = t — p. 
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FIG. 4. Relationship of r, p and ^, p coordinates. The dark bars represent the peaks of the X oscillations; 
the dashed line shows the location of the transition edge. In both coordinate systems the region well above 
the transition edge is the "quiescent" region, and the region well below the edge is the "oscillatory" region. 



Figure 4 shows schematically the location of the quiescent and oscillatory regions and of the 
transition edge, in terms of the r, p and r, ji coordinates. In terms of these last coordinates, the 
Einstein-scalar equations are, 

-aY + {aY)^p~{aY)^^, (15) 
aX + {aX)^p-{aX)^^, (16) 
X^+y^ (17) 

2a^{X'^ + Y'^) . (18) 

As an example of the usefulness of the p coordinate, we now go back to the behavior of a in (|lj) 
and consider the nature of T{t). At any r the change in the "source" a^{X^ + Y'^) becomes sizeable 
at the transition edge p k, Q. This tells us that at different times the characteristic rise of a as a 
function of p will always occur around /i w 0. For this to be true, and for a to have the form e^^ at 
any one r, requires that the dependence of a on r and p have the form 

a-e"''^ (19) 

in the oscillatory region, i.e., for p <^ —\. This approximation gives a a self-similar form, and is 
too crude, since the exact location at which X, Y start to be large is not really p = Q but (since the 
solutions only have a discrete self similarity) oscillates in r. A better approximation which takes 
this oscillation into account, is 

a = A(T)e^'"^ , (20) 
where A{t) is a bounded, oscillatory function. 
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III. APPROXIMATE SOLUTIONS 



A. Oscillatory region /i << — 1 

We make the annate (motivated by the resuhs in Figs. 1 and 2) that the /x-derivative terms on 
the right in equations (^5|) and (|l^) are to be ignored in comparison with the e^^9^ derivatives on 
the left. We also introduce the variable 

w = t/r = -e^ (21) 

in terms of which the equations take the form 

(aX)^„ = -aY + {aY)^p (22) 
{aY)^u ^aX + {aX)^p . (23) 

If we treat a and a as constant, the equations admit the solution 

Y ^ Aii\n[uj{p - au)] (24) 

X = A sin[cj(p - au) ~ 5] (25) 

with, 

l + ^^a^ (26) 

cos(5 = -^ (27) 

sin^ = — . (28) 
uja 

If one fits this solution to the Choptuik's data in Figs. 1 and 2, one finds that 5 ~ 150 degrees. 
The ratio of side lengths of the phase space ellipse is given by tan (5/2), which yields a value of 
approximately 4, agreeing with the data due to Choptuik exhibited in Fig. 3. 

Since we have assumed that a and a are constant, we have in effect been considering the spacetime 
to be a flat background. This suggests that we should be able to find similar behavior in a solution 
to the scalar field wave equation in Minkowski spacetime. In this case </) can always be expressed, 
in terms of arbitrary functions /, g, as 

r 

In terms of / and g, Choptuik's phase space variables X, Y are 

Y = f\t-r)+g'{t + r) (30) 

X = --{}{t-T)+g{t + r))- nt-r)+g\t + r) . (31) 
r 

To match a solution of this form to that of (|2^), ( p5| ) we take / and g to have forms chosen to give 
good agreement with ( p^ , (|5|) for /i <C — 1 ; 

\ -\- k Tj 

f(.v) = ^^r~ — 9 In 0V) - ^ cos {uj In f3r])] (32) 

2 1 + CiJ^ 

g{ri) — A — ^ [sm {lu \n f3ri) — UJ cos {uj In I3rj)] . (33) 

Here A and /? are adjustable constants and — 1 + l/cu'^. The agreement of this solution with the 
numerical results of Choptuik is most impressive in a phase-space plot, as in Fig. 5, where we have 
used the values A = 0.24, a; = 1.85, A; = 1.137 in the flat spacetime solution. (The value of /3 does 
not infiuence the appearance of the solution in this plot.) 
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FIG. 5. A phase-space plot comparing Choptuik's numerical results with a flat spacetime solution at 
T = —7.87488. The flat spacetime solution corresponds to parameters uj = 1.85 and A — 0.24. 

The light dotted curve is Choptuik's result for t = —7.87488 and the thick curve is the flat 
spacetime solution. The thick curve, of course, lacks the vertical spike representing the transition 
edge. For the flat spacetime solution there is no transition front. The solution in (]3^),(^3|) is a valid 
solution except ati-|-r = Oort — r = (and is shown in Fig. 5 only for t — r < and t + r > 0). 
It can agree with the Choptuik numerical results only in the region in which the approximations 
(p^) , (|2^) apply, i.e., it can agree only in the oscillatory region. 

The flat spacetime fit of eqs. (|3^) , (p3|) , is plotted and compared to the numerical results of Chop- 
tuik in Figs. 6a, b. The agreement is remarkable in the oscillatory region, except at large p (where 
Choptuik's results deviate from the ideal critical solution). The flat spacetime solution stops being 
reasonably accurate, of course, near the transition front, and deviates wildly from Choptuik's results 
in the quiescent region. 

Despite the impressive agreement (in the appropriate region) of the simple solution above with 
the numerical results, the solution is less than ideal. For one thing it does not explain the gradual 
decrease in the amplitude of the oscillations seen in the numerical results. More worrisome is 
the question whether it is justified to treat a as constant when Fig. 2 shows it to be increasing 
exponentially. An improved approximation scheme which deals with such issues is possible, but is 
not central to our main point of this paper. 
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FIG. 6. The solutions X{p), Y{p) found by Choptuik for r — —7.87488 are compared to the flat spacetime 
solutions corresponding to parameters u — 1.85, A — 0.24, l3 = 0.35 



B. Quiescent region /i >> 1 

In this region we drop the terms involving e~''9p on the left of and (p^), and keep the 9^ 
terms. We make the convenient definitions x = aX and y = aY so that the equations for X and Y 
(with the e~^d^ terms dropped) read, 

J/,p - y,^ = y (34) 
x.p X — X . (35) 

The general solution to these equations can immediately be written as 

y ^ e'^-f-+p)/^F{T) X ^ pJ^-^^/^GCr) , (36) 

in which F and G are arbitrary functions. At /i = these solutions become 

y = eP'^F(p) X = e-p'^G{p) . (37) 
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FIG. 7. The solutions X{p), Y{p) found by Choptuik for r — —7.87488 are compared to the quiescent 
region flat spacetime solutions. Parameters A = 0.324, 5 — 0.372, have been used in the flat spacetime 
solution given in the text. 

But at /i = 0, the solution is bounded a,s p ^ —oo, so F{^) and G{^) must respectively have the 
form e^''/^F(^), e^/^G'(^), where F(^), G'(^) are bounded as ^ ^ — oo. In the quiescent region, then, 
X and y have the form 



y 

X 



Fir) 
G{t) 



(38) 
(39) 



At constant r, in the quiescent region, the scalar field variables vanish as p — » — oo. This is satisfied 
by the above form for y, but for x not to diverge, G must be zero. This "first approximation" to 
X and y shows why X/Y — > as p — > — oo, and why the phase space plots have vertical spikes 
corresponding to the p — > — oo end of the plot. We can find a better solution for x by substituting 
the above solution for y into the left hand side of (|l6|), and by treating a and a as constants. This 
leads to 



H{fi + p) , 



(40) 



where H is given by 
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a 

3^ 



(41) 



From the results of Choptuik, such as those shown m Figs. 1 and 2, it is clear that the solutions 
in the quiescent region oscillate in t, with a period of 27t/lu. The function F, then, must have this 
periodic property. 

In the above, we have considered a and a to be constant, so the solutions in ( |3^ ) and (^0|), must 
be flat spacetime solutions. To see this explicitly, we use the flat spacetime forms of the solution in 
( pO| ) and (pl|). Since Y is finite as r ^ we must have /(yy) = —g{ri). In the limit of small r this 
gives us 



Y = -2rf"{t) + 0{i 



(42) 



By comparison with ( p8| ) we have that f"{t) — —F{T)/2t, where F is a periodic function with period 
2n/uj. We illustrate the agreement of the Choptuik results with such a flat spacetime solution, by 
choosing 



/(ry) = In 77 sin {ujlnri + 6). 



(43) 



Figure 7 shows the agreement, in the quiescent region, between the Choptuik results and the flat 
spacetime results of (jl|) , (|l6l) , (||) , at r = -7.87488. In (|3|) the parameter values A = 0.324,(5 = 
0.372, have been used. 



IV. DISCUSSION 



In the previous section we introduced two approximate solutions that describe, with reasonable 
accuracy, the behavior of the Choptuik spacetime and we showed that these solutions were equivalent 
to flat spacetime solutions. Neither of these flat spacetime solutions, by itself, provides an acceptable 
solution in the full range of the spacetime manifold in which the Choptuik solution is found. From 
this point of view, the essential role of general relativity, and the nonlinearities it introduces into 
the equations, is to provide a match of the two solutions across the transition edge. We give here 
a very rough version of such a match. This is done to illustrate orders of magnitudes, but it is 
amusing (and coincidental) that this approximate match ends up giving a fairly accurate estimate 
of the period of the Choptuik echoing. 

We know from the flat spacetime approximation in the echoing region that, 

1 + . (44) 

We now adopt the viewpoint that to is fixed by the way in which a increases across the transition 
edge. The requirement of elementary flatness fixes a to be unity at at r = 0. Since X and Y are 
small in the quiescent region, we see from (^, here rewritten as 

(-).p + - = x2 + r2^ (45) 

r r 

that m/r must remain very small in the quiescent region, since X and Y are small. It follows that a 
must be very close to unity in the quiescent region, and this is verified by Choptuik's results in Fig. 
2. Equation ([45|) also allows us to make a crude estimate of how a increases across the transition 
edge. From (pq) the discussion in the previous section we can approximate X Ri and (with a taken 
to be approximately unity) Y « e^^i^(r). This approximation in ( ^5[ ) gives us 

{-),, + -=e'''e-^^F\r). (46) 
r r 

The initial condition on the solution is that m/r = at the origin (p — —00), so that the solution 
everywhere in the quiescent region is approximately 

- = \F'ir)e~^^ . (47) 
r 6 
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We now assume that this solution gives a roughly correct answer at the transition edge, fj, = 0. Since 
F is an oscillating function, m/r also oscillates (see Fig. 2), and a representative value of m/r is 
given by taking an average: 

<!!^>;^i<F2> . (48) 
r S 

We now argue that F must be of order unity. The transition edge is defined as the location at which 
X and Y grow to order unity. If this is to be at /i = then, from (^sj) it is clear that F must be 
^ 1. If, in fact, we take to be a sinusoidal oscillation of amplitude unity, then < F^ >« 1/2, and 
we have from (Ea) that < m/r >« 1/6. 



As a consequence, 



which together with (MJ) implies. 



9 2m 4 , , 

a2 ^ 1+ < _ >^ - , (49) 



r 



- Va, A ^ - 3.6 (50) 

v3 

Of course, we could equally well have concluded, for example, that < m/r >« 0.1 and A ~ 2.8. Our 
argument can only establish that A is of order unity. But this is enough to support our suggestion 
that it is the matching provided by the nonlinearities that picks out the form and the period of the 
Choptuik solution as the only one that can bridge the two approximately flat-spacetime solutions. 
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